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We extend the method of time delay proposed by Eisenbud and Wigner, 
l/") ' to search for unstable states formed by n mesons and the 3 He nucleus. 

s ! , Using few body equations to describe 7)- 3 He elastic scattering, we predict 

resonances and unstable bound states within different models of the TjN 
interaction. The n 3 He states predicted within this novel approach are in 
agreement with the recent claim of the evidence of ry-mesic 3 He made by 
the TAPS collaboration. 
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A claim of positive evidence for the existence of the long sought after 77-mesic 
nuclei was made last year by the TAPS collaboration pQ on the basis of their 
investigations on the photoproduction reaction, 7 3 He — ► r/X. This indeed comes 
"q ■ as a big step in the understanding of many nuclear and particle physics related 

^ I problems of fundamental interest. We list here the connection of 77-physics to 

a few of these important issues: (i) the 77-N interaction is dominated by the 
coupling to the A?"* (1535) resonance and hence, studies of the 77 meson in the 
nuclear medium could yield information on the r]NN* coupling constant, (ii) 
The isospin and charge symmetry breaking (CSB) in quantum chromodynamics 
(QCD) arises due to the difference in the up and down quark masses due to their 
electromagnetic interaction. The observation of the ir-rj-r)' mixing provides one 
of the best possibilities to study directly the CSB effects [2]. (hi) The relative 
cross sections for the reactions np, rj'p — > rjp, rj'p and 7r~ p — > 7771, rfn, provide 
a sensitive test for the presence of a strange-antistrange (ss) component in the 
nucleon's wave function Obviously then, this evidence of an ?7-mesic 3 He 
increased the interest in this field which started about a decade ago |S] . We 
therefore found it timely to search for the 77-mesic resonances using the method 
of time delay proposed by Eisenbud and Wigner [S] about half a century ago, 
but revived only recently to successfully characterize the meson and baryon 
resonances |7||S]. From our experience in [7], we found that with this method 
(which is mentioned in text books and literature as a necessary condition for the 
existence of a resonance) we not only confirmed the existing resonances, but also 
found new ones. Our finding of the exotic pentaquark resonance around 1540 
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MeV from K + N data [5] was indeed confirmed by several recent experiments 
(listed in 

The time delay in rj 3 He elastic scattering is evaluated using few-body equa- 
tions for the T\ 3 He system involving inputs from different models of the ele- 
mentary r/N interaction. In an earlier work UJ, the near threshold data on the 
pd — > 3 He rj reaction were well reproduced using the same few body equations 
to describe the 7? 3 He final state interaction. In what follows, the method of 
time delay is briefly introduced, followed by an example which demonstrates 
the usefulness of this concept even at negative energies. 

Back in the early fifties, Eisenbud and Wigner related the energy derivative 
of the phase shift to the time delay in scattering as, At(E) = 2hd5/dE. A 
more useful definition, namely, the time delay matrix was given by Eisenbud 
and elaborated by Smith ^Uj. An element of this matrix, Atij, which is the 
time delay in the emergence of a particle in the j th channel after being injected 
in the i th channel is given by, 



where Sij is an element of the corresponding S-matrix. In an eigenphase formu- 
lation of the S'-matrix, one can easily see that the time delay as defined above 
in is the energy derivative of the phase shift. Alternatively, writing the S- 
matrix in terms of the T-matrix as, S = 1 + 2 i T, one can evaluate time delay 
in terms of the T-matrix. The time delay in elastic scattering, i.e. At a, is given 
in terms of the T-matrix as, 

S* Su At u (E) = 2h 5Re + 2 ReT« 3m^j-2 3mT„ ; Sfte 

'(2) 

where T is the complex T-matrix such that, Tkj — tfteTkj + i^smTkj. The 
above relations were put to a test in |7| to characterize the hadron resonances 
occurring in meson-nucleon and meson-meson elastic scattering. The energy 
distribution of the time delay evaluated in these works, nicely displayed the 
known N and A excitations and meson resonances like the p, scalars (/o) and 
strange if*'s found in Kit scattering, in addition to confirming some old claims 
of exotic states. Evidence for the low-lying exotic pentaquark was also found in 
addition to several other Z*'s from the time delay in K + N scattering. 

We shall now demonstrate that the time delay concept which has been so 
useful in characterizing resonances, can also be extended to the search of bound 
as well as unstable bound states. The latter signify states with a negative 
binding energy but having a finite lifetime. These are sometimes called as 
"quasibound" states in literature. We consider the example of an S'-matrix for 
the n — p system written for the case of a square well potential which has the 
right parameters to produce the binding energy of the deuteron. The S matrix 
in this CclSC, 0jS cL function of / is given as, 

ahf> \a)M0) - /3^ (1) (")i;(/3) 
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Figure 1: The theoretical delay time in np elastic scattering as a function of 
the energy E = \fs — m n — m p , where \fs is the total energy available in the 
np centre of mass system. The sharp spike at E = —2.224 McV corresponding 
to the deuteron binding energy indicates the infinite time delay due to the 
formation of this np bound state. The shaded distribution corresponds to the 
time delay in the case of a fictitious unstable bound state. 



where ji , h\ and h\ are the spherical Bessel and Hankel functions of the first 
and second kind respectively, a = kR and (3 = (a 2 — 2[iU R 2 /h 2 ) 1 / 2 where the 
potential U is given by — U = V + iW , with U(r) = U 8(R — r) and R the width 
of the potential well. A similar square well potential was used by Morimatsu and 
Yazaki while locating the "unstable bound states" of E-hypernuclei as second 
quadrant poles and by J. Fraxedas and J. Sesma using the time delay method 

EH 

We evaluate the time delay in n — p scattering at negative energies using 
the above ^-matrix with a — ikR (hence E = — fc 2 /2/i) and I = 0. If we plug 
in the appropriate parameters for an n — p square well potential, V = 34.6 
MeV and R — 2.07 fm, the time delay plot (as in Fig. 1) shows a sharp spike 
exactly at the binding energy of the deuteron (E = —2.224 MeV). If we add 
a small imaginary part to the potential, then this state develops into a Breit- 
Wigner kind of distribution centered around the binding energy of the deuteron 
(a fictitious "unstable bound state" of the n — p system around E = —2.224 
MeV). We can thus see that time delay is well-defined at negative energies 
and does reveal the bound and "unstable bound states" with positive peaks at 
negative energies. At this point, we may add that in contrast to the infinite 
time delay that occurs in the case of bound states, the delay time for virtual 
states (also called "antibound" states) is — oo. Below, we describe our approach 
to construct the complex T-matrix for r\ 3 He elastic scattering which is required 
for the calculation of time delay. 



3 



The 77 3 He transition matrix is evaluated using few body equations solved 
within a Finite Rank Approximation (FRA) approach, which means that the 
3 He nucleus in fy 3 He elastic scattering remains in its ground state, in the in- 
termediate state. Since the 77-mesic bound states and resonances are basically 
low energy phenomena, it seems justified to use the FRA for calculations of the 
present work. Since at low energies, the 77.ZV interaction is dominated by the Sn 
resonance, we restrict ourselves to the study of the s-wave unstable states. The 
?7 3 He i-matrix in this approach is written as |12l 

tr, A (k',k;z) = < 9; ip \t°(z)\k; > + (4) 
dk" < k';^ \t°(z)\k ! ';^ > ^ 7/ -> 
(2-) 3 (z-tgXz-e-tg) tvMk ' *''*>' 

where z = E — \e\ + iO. E is the energy associated with r\A relative motion, e is 
the binding energy of the nucleus, ipo is the nuclear wave function and /1 is the 
reduced mass of the r]A system. The matrix elements for t° are given as, 

< k'; ^ Q \t°(z)\k; V-o >= J dr\Mr)\ 2 t°(k'J;r;z) (5) 

where, 

A 

t°(k',k;r;z)=J2ti(k',k;n;z). (6) 

i=l 

i° is the t-matrix for the scattering of the fy-meson from the i th nucleon in the 
nucleus, with the rescattering from the other (A-l) nucleons included. It is 
given as, 

/ 1~ */ j* — * \ tjN / ' iT/ r* — * \ f dk ^/c . Jx . 1**1 . z^j ^ / 7"// J* — * \ 

t$ (re , re ; Ti ; z) = t i (k , A; ; r-j ; z)+ / . . 3 > * > r i ! z ) ■ 

J [ K) Z 2/T jjti 

(7) 

The t-matrix for elementary ry-nucleon scattering, t^ N , is written in terms of 
the two body rjN matrix i^N^r/Af as, 

tV (re', k;fl;z) = t vN ^ vN (k', k ; z) exp[i( k - k' ) ■ fl) . (8) 

Since there exists a lot of uncertainty in the knowledge of the 77-nucleon interac- 
tion, we use three different prescriptions 1141 1TB] of the coupled channel 77-N 
t-matrix, t v n —>n N '■> leading to different values of the rjN scattering length. 

Before proceeding further, we briefly list the general features of the present 
approach, the approximations made and their validity in the context of the 
present calculations: 

(i) The main idea of the FRA consists of separating the motion of the projectile 
and the internal motion of the nucleons inside the nucleus. The total Hamil- 
tonian is split accordingly and the approximation in solving these equations as 
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mentioned above is to restrict the spectral decomposition of the nuclear Hamil- 
tonian to the ground state. This, however, limits the applicability of the present 
approach to energies below the break up threshold of the 3 He nucleus. 

(ii) The operator t° describes the scattering of the ij meson from nucleons fixed 
in their space position within the nucleus. However, it differs from the usual 
fixed center t-matrices as t° is taken off the energy shell and involves the motion 
of the r\ meson with respect to the center of mass of the target. The present 
scheme should not be confused with an optical potential approach which in- 
volves an impulse approximation and the omission of higher order rescattering 
terms. Besides, the i-matrix is evaluated at z = E — \e\ + iO in contrast to the 
usual fixed scatterer approximations which do not involve the binding energy 
term. 

(iii) the dominance of the N* — Sn resonance which decays with a similar 
branching ratio to both the ttN and ijN channels implies that the transitions 
77 N — ► irN must be taken into account. In the present work we use coupled 
channel t-matrices for the elementary ^-nucleon interaction which include fully 
the effect of the ttN and rjN channel on each other. The complex self energy 
terms in the r]N — > r]N t-matrix take care of the intermediate off shell as well 
as on shell ttN loops. We neglect however, transitions involving the production 
of an rj meson on a second nucleon by a pion produced on the first nucleon. For 
example, in the present three nucleon case, we neglect processes of the type, 
■qNi — ► irNi followed by 11N2 — > r/N2 or ttN% — ► ttN2 followed by itN% — ► 77^2 
etc. These are in principle more difficult to calculate. Though this neglect 
cannot be justified a priori, there is reason enough to believe that these contri- 
butions may not be significantly large as argued in [16] . 

(iv) The 3 He nuclear wave function required in the calculation of the T-matrix is 
taken to be of the Gaussian form. The use of this wave function for the present 
work seems adequate in the light of the calculations performed in 13^ where 
searches for eta-mesic states using a more sophisticated wave function for 3 He 
resulted in almost the same locations of the T-matrix poles as found using the 
Gaussian one. 

The T-matrix for r\A elastic scattering, t v A, is related to the S'-matrix as, 

S = l- f —t nA , (9) 

TT 

where k is the momentum in the rjA centre of mass system and hence, the 
dimensionless T-matrix required in the evaluation of time delay is given as, 
T = — (fi kj 2 7i") t v A- and used to evaluate the time delay given in © for the re- 
action 7? 3 He — > 77 3 He. We evaluate Atu(E) at both positive and negative values 
of E in order to search for resonances and unstable bound states. As explained 
above with the example of the deuteron, time delay at negative energies is a 
valid concept and is useful in identifying the bound and unstable bound states. 
However, in contrast to positive energy resonances where one can "measure" 
the time delay from phase shifts in a scattering experiment, the time delay at 
negative energies is not measurable in an elastic scattering experiment. At neg- 
ative energies however, it manifests in an off-shell elastic scattering, where the 
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transition matrix t^A is calculated at purely imaginary momenta (k = k' — ► ik) 
corresponding to E = —k 2 /2/j, with /i being the reduced mass of the rj A system. 
Physically, this would correspond to the binding of subthreshold produced 77's 
in nuclear reactions through off shell scattering. 
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Figure 2: Time delay in the elastic scattering reaction, 77 3 He — * ?? 3 He, as a 
function of energy, E = E„3jj e — m n — TO3# e , where E v 3H e is the total energy 
available in the ?y 3 He centre of mass. The vertical axis scale has been broken 
to display the structure of the time delay plot clearly. 



In Fig. 2, we show the ?7 3 He time delay plots with two different inputs for 
the r)N interaction In ^3] a coupled channel t-matrix including the 

7rN and ?yN channels with the Sn - 77N interaction playing a dominant role was 
constructed. It consisted of the meson - N* vertices and the N* propagator as 
given below: 

t v N ^ N ( k', *; z) = r N , (z) J^Ly (10) 

with, t n ,(z) = (z - M - T,^(z) - E v (z) + ie)~\ where T, a (z) (a = tt, 77) 
are the self energy contributions from the irN and r)N loops. The parameters 
were fitted such that they reproduced scattering lengths of a vN = (0.75, 0.27) 
fm and a^jv = (0.88, 0.41) fm, which are in agreement with values obtained 
from modern analyses of recent data. We chose the parameter set with g N „ — 
2.13, P = 13 fm" 1 and M = 1656 MeV, which gave rise to a^N = (0.88, 0.41) fm. 
In 0], the 7rN, 77N and irA (ttttN) channels were treated in a coupled channel 
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formalism (so that an additional self-energy term appears in the propagator 
t n ,(z)). The parameters of this model are, g N , — 0.616, (3 = 2.36 fm _1 and 
Mo = 1608.1 MeV. These parameters were obtained from a fit to the P33 and 
S11 phase shifts in ttN elastic scattering and the differential irN — > r/N cross 
sections predicted by the model were in good agreement with data. The input 
parameters of this model, however, give rise to a much smaller scattering length, 
namely, a V N = (0.28, 0.19) fm. It should be noted that though the two models 
predict somewhat different total ttN — >■ rjN cross sections, they both agree with 
data since the error bars on the data are large. 

In both cases we find a peak centered slightly below zero energy. In the case 
of the larger rjN scattering length ^3], we find an additional bump located at a 
positive energy of 0.5 MeV and overlapping the peak very close to threshold. In 
the case of the smaller a V N (lower half), there is a hint of a broad hump around 
—5 MeV. These findings seem to be in agreement with experiment where a 
binding energy of — 4.4 ± 4.2 MeV for the ry-mesic 3 He and a resonance like 
structure just above production threshold was reported. 

Before proceeding further, we would like to mention that in the case of s- 
wave scattering near the elastic threshold, i.e. at zero energy, one needs to be 
careful in drawing conclusions from time delay plots. Using S = exp(2iS) and 
comparing it with 

8 = i ln(l - = I ln(l + 2ikf) (11) 

Zl 7T Zl 

where / is the scattering amplitude. For small k, 8 ~ kf and the behaviour 
of d8/dE (the real part of which is essentially the time delay) is determined by 
the simple pole at k = (or E v 3 He = E t hreshoid) and the energy dependence of 
the scattering amplitude /. In the absence of a resonance, as k — > 0, 8 = ka, 
where a = an + iaj is the complex scattering length. For positive energies, 
$ie8 = kaii, whereas for energies below zero, k — > ik and 5ie<5 = —kaj. In such 
a situation, $te(d8 / dE) exhibits a sharp peak at E v 3jj e = Ethreshoid, the sign of 
which is determined by the sign of the scattering length. On the other hand, if 
the scattering amplitude has a resonant behaviour near threshold, one would see 
a superposition of the two behaviours. For a state far from threshold, the rise in 
time delay near threshold is extremely sharp and the shape of the state remains 
completely unaffected. A state reasonably close to threshold is distorted in shape 
and one very close manifests simply by broadening the threshold singularity. In 
Fig. 2, we display the phase shifts and their threshold behaviour (dashed lines) 
which depend on the rj— 3 He scattering lengths. These phase shifts do display 
the characteristic resonant behaviour apart from their agreement with 8 = ka 
close to threshold. 

The formalism of the present work forces the three nucleon system to be in 
the bound state always and hence searching for unstable states beyond the 3 He 
break up threshold of 7.7 MeV as well as below -7.7 MeV is beyond the reach 
of the present approach. This limitation taken along with the fact that virtual 
target excitations could be important below the break up threshold add some 
uncertainty to the bump around —5 MeV. However, this bump is quite broad 
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and spreads over to the positive energy side near threshold where the FRA is 
reliable. Besides, we note that the r]N scattering length in this case is very small 
(a n N — (0.28,0.19) fm) and in a comparison of the Alt-Grassberger-Sandhas 
(AGS) equations (which include the target excitations) with FRA, the authors 
in J7| found that the FRA works reasonably well for real parts of a V N < 0.5 
fm. The above peak is also expected to be distorted by the singularity in the 
threshold region as mentioned above. Thus, one cannot be sure about the exact 
peak position of this state, but the plot does seem to indicate the existence of 
one broad state in this region within the eta-nucleon model of Bhalerao and 
Liu. 

Another interesting feature in these time delay plots is the occurrence of 
negative time delay around 7.7 MeV which is essentially the binding energy of 
the 3 He nucleus and is also an input to our calculation (see Eq. 0}. It can be 
understood by noting the connection of the density of states with the energy 
derivative of the phase shift which also defines time delay. According to the 
Beth-Uhlenbeck formula |18|. the difference in the density of states with and 
without interaction is given by the energy derivative of the phase shift. This 
means that when the density of states due to interaction is less than that without 
interaction, the energy derivative of the phase shift would be negative. The 
reduction in the density of states due to interaction in our context corresponds 
to a loss of flux from the elastic channel due to the opening of an inelastic 
channel (more examples in |19p. 

In Fig. 3, we show the time delay plots and corresponding phase shifts using 
a simple separable form of the off-shell t-matrix ^H] . Though this t- matrix does 
not contain the self-energy contributions as in 0] ^] > it has the advantage that 
it is parametrized with the latest Crystal Ball data on n N — ► i]N [23] and also 
includes the 7-/V — > r/N data In Fig. 3a, we observe subthreshold peaks 

for two sets of t-matrix solutions (Set A and B from |15|1. out of which the 
solid line indicates their best solution and dashed one a similar solution without 
pion beam momentum correction as explained in |15) . In both these cases, 
however, we observe negative time delay above threshold which could be due to 
a repulsive interaction ^J^l- I n Fig- 3b, one sees a broad peak at threshold, 
very similar to those in Fig. 2. This result corresponds to Set D in ^j] which is 
an unconventional solution obtained when the recent data from |20| is omitted. 
Finally, Fig. 3c displays two distinct peaks using the Set E solution which 
omits the data from Though we have listed only the scattering lengths 

corresponding to the various solutions, it should be noted that the results are 
also sensitive to the effective range parameter ro given in Table III of |15| . 

Before concluding, we mention an analysis 22 where, using optical poten- 
tials, the authors investigated the relation between the ry- 3 He binding energy, 
width and the complex scattering length, a, ( _3 ffe . The strength of the potential 
was varied to check for which scattering lengths one can find bound states. They 
found that if an 77-mesic 3 He state exists, its binding energy (width) should not 
exceed 5 MeV (10 MeV). It is interesting to note that though the approach of 
the present work which involves the indentification of the quasi bound state from 
positive peaks in time delay is entirely different from that in [22], the possible 
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Figure 3: Same as Fig. 2, but with a different model of the r]N interaction. 
Scattering lengths are given in fm. 



77-mesic states of the present work lie within this limit. The constraints on the 
energy and width in (22| were found following a systematic analysis of the ry- 3 He 
scattering length ^23: . The values of a v _3 He obtained in the present work are 
outputs of the few body calculations for a given rjN interaction that enters the 
77-nucleus i-matrix as an input. They are listed in Figs 2 and 3. We refer the 
reader to [231 an d the references therein for a detailed discussion on the issue of 
the 77-nucleus scattering lengths. 

In conclusion, we have investigated the possible existence of 77-mesic 3 He us- 
ing the physical concept of time delay and within the limitations of the present 
approach, find evidence for the existence of such states near threshold using 
two off-shell isobar models (|J and jTl]) of the r/N interaction and another 
if -matrix semiphenomenological analysis of relevant data With the avail- 
ability of more accurate data on 77-mesic nuclei, the present results will be useful 
in obtaining a better knowledge of the rjN interaction. 

References 

[1] M. Pfciffcr et at, Phys. Rev. Lett. 92, 252001 (2004). 

[2] A. Magiera and H. Machner, Nucl. Phys. A674, 515 (2000). 

[3] C. B. Dover and P. M. Fishbane, Phys. Rev. Lett. 64, 3115 (1990). 



9 



[4] R. S. Bhalerao and L. C. Liu, Phys. Rev. Lett. 54 (1985) 865. 

[5] Q. Haider, L. Liu, Phys. Lett. B 172, 257 (1986); ibid 174 (1986) 465(E). 

[6] E. P. Wigner, Phys. Rev. 98, 145 (1955); E. P. Wigner and L. Eisenbud 
Phys. Rev. 72, 29 (1947); L. Eisenbud, dissertation, Princeton, June 1948 
(unpublished) . 

[7] N. G. Kelkar, M. Nowakowski, K. P. Khemchandani and S. R. Jain, Nucl. 
Phys. A730, 121 (2004), lh^ph/ 0208197. N. G. Kelkar, M. Nowakowski 
and K. P. Khemchandani, Nucl. Phys. A724, 357 (2003). 

[8] N. G. Kelkar, M. Nowakows ki and K. P. Khem chandani, J. Phys. G: Nucl. 
Part. Phys. 29, 1001 (2003),|hep-ph/0307134| ibid, Mod. Phys. Lett. A19, 
2001 (2004), |nuchth/0405008| ' 

[9] K. P. Khemchandani, N. G. Kelkar and B. K. Jain, Nucl. Phys. A708, 312 
(2002), IrmcTth/Ol 12065| 

[10] F. T. Smith, Phys. Rev. 118, 349 (1960). 

[11] O. Morimatsu and K. Yazaki, Nucl. Phys. A435, 727 (1985); J. Fraxedas 
and J. Sesma, Phys. Rev. C 37, 2016 (1988). 

[12] V. B. Belyaev, Lectures on the theory of FEW BODY SYSTEMS, Springer 
series in Nuclear and Particle Physics, 1990. 

[13] S. A. Rakityansky et al, Phys. Rev. C 53, R2043 (1996); S. A. Rakityansky, 
S. A. Sofianos, V. B. Belyaev, W. Sandhas, Few-Body Syst. Suppl. 9, 227 
(1995). 

[14] A. Fix and H. Arenhovel, Eur. Phys. J A 9, 119 (2000); ibid, Nucl. Phys. 
A 697, 277 (2002). 

[15] A. M. Green and S. Wycech, Phys. Rev. C 71, 014001 (2005). 

[16] A. Fix and H. Arenhovel, Phys. Rev. C 68, 044002 (2003). 

[17] N. V. Shevchenko et al, Phys. Rev. C 58, R3055 (1998). 

[18] E. Beth and G. E. Uhlenbeck, Physics 4, 915 (1937); K. Huang, Statistical 
Mechanics (Wiley, New York, 1963). 

[19] N. G. Kelkar, J. Phys. G: Nucl. Part. Phys. 29, LI (2003), piep^ph/0205188| 

[20] N. G. Kozlenko et al, Gatchina Preprint 2542, 2003 (unpublished). 

[21] F. Renard et at, Phys. Lett. B 528, 215 (2002). 

[22] A. Sibirtsev et al, Phys. Rev. C 70, 047001 (2004). 

[23] A. Sibirts ev et al, Eur. Phys. J A 22, 495 (2004); e-Print Archive: 
|nucl-th/0310079| 



10 



